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I  INTRODUCTION 


This  report  describes  the  results  of  an  investigation  of  crater 
formation  by  nuclear  surface  bursts  done  during  the  1978  and  1979  JASON 
Summer  Studies.  Early  time  processes  (radiation  and  hydrodynamics)  and 
scaling  laws  are  emphasized. 


1 . 1  General  Remarks 

Before  discussing  in  detail  the  contents  of  this  report,  let  us 
first  make  clear  the  spirit  in  which  this  work  is  presented.  We  have  not 
attempted  a  complete  theoretical  or  experimental  review  of  the  nuclear 
cratering  problem;  nor  have  we  undertaken  numerical  computer 
calculations.  Rather  our  report  is  intended  to  be  a  simple  presentation  of 
the  important  physical  processes  involved  in  nuclear  cratering  along  with 
their  spatial  and  temporal  scales.  We  hope  that  this  will  be  helpful  in 
interpreting  the  results  of  the  large  scale,  state-of-the-art  numerical 
calculations  and  also  in  resolving  some  of  the  disagreements  and 
controversies  which  currently  exist.  We  recognize,  of  course,  that  the 
definitive  theoretical  work  on  this  subject  will  most  likely  come  from 
large  computer  studies  which  model  in  detail  the  complex  phenomena 
involved.  We  recognize  too  that  the  existing  database  for  craters  formed 
by  surface  bursts  of  nuclear  weapons  is  very  limited  and  subject  to  a 
variety  of  interpretations.  Much  of  what  we  say  is  probably  already 
contained  in  the  notebooks  of  workers  in  the  field.  We  hope  that  this 
report  will  encourage  a  wider  circulation  of  simple  formulas  and  order  of 
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magnitude  calculations,  so  that  a  more  complete  cratering  "primer"  can  be 
built  up  and  a  convergence  of  opinion  achieved  on  the  important  mechanisms 
of  this  complex  problem. 

Finally,  let  us  mention  our  point  of  view  concerning  the  overall 
role  of  nuclear  cratering.  For  most  offensive  and  defensive  purposes  the 
dominant  mechanism  of  energy  coupling  into  the  ground  in  surface  and  near 
surface  bursts  is  the  airblast  induced  ground  shock,  not  the  direct  ground 
shock  which  is  co-produced  with  the  crater.  The  distance  beyond  which  the 
relatively  well  understood  airblast  phenomena  dominate  the  more  poorly 
understood  direct  phenomena  is  on  the  order  of  a  few  crater  radii.  Thus, 
trustworthy  models  for  crater  characteristics  and  related  phenomena  are 
required  for  only  a  very  restricted  set  of  applications  typically  involving 
superhard  or  deeply  buried  targets.  We  have  not  studied  the  nature  and 
merits  of  such  applications  but  rather  have  accepted  them  as  given. 

1.2  Contents 

We  now  list  the  contents  in  this  report  and  then  summarize  our 
conclusions  and  recommendations.  In  section  II  we  give  a  brief 
introduction  to  the  physical  phenomena  which  are  involved  in  crater 
formation  in  nuclear  surface  bursts  (i.e.,  radiation,  hydrodynamics, 
elasto-plastic  flow,  etc.).  In  section  III  we  present  a  discussion  of 
scaling  laws  which  have  been  applied  to  cratering.  We  begin  with  a  review 
of  simple  dimensional  analysis  and  then  discuss  some  recent  findings  from 
centrifuge  cratering  studies.  We  then  review  the  phenomenological  scaling 
laws  which  have  been  abstracted  from  available  data  and  make  a  few 
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comparisons  to  the  centrifuge  findings.  Lastly  we  summarize  comparisons 
between  presently  available  computer  calculations  and  nuclear  craters 
formed  by  surface  bursts. 

In  section  IV,  we  present  a  crude  calculation  of  the  radiation 
coupling  to  the  ground  and  the  initial  impulse  thereby  generated.  This 
topic  is  a  particularly  critical  one,  since  most  of  the  energy  release  from 
modern  nuclear  weapons  comes  in  the  form  of  a  short  burst  of  kilovolt 
X-radiation.  The  controversial  subject  of  fireball  coupling  is  discussed 
here . 

In  section  V,  we  discuss  the  hydrodynamic  era  in  which  the  initial 
impulse  produced  by  the  radiation  induced  blow-off  is  amplified  by  further 
ejection  of  shock  processed  matter.  Here,  as  in  the  preceding  section,  the 

| 

presentation  is  a  "back  of  the  envelope"  one  aimed  at  order  of  magnitude 
estimates.  The  subsequent  elasto-plastic  flow  is  addressed  in  only  a  crude 
manner. 

A  final  section  (VI)  presents  our  conclusions  and 
recommendations.  These  are  based  in  part  upon  the  order  of  magnitude 
calculations  presented  in  the  earlier  sections  and  in  part  upon  the  written 
and  oral  presentations  we  received  prior  to  and  during  our  studies. 


II  PHYSICAL  PHENOMENA 


We  present  here  a  brief  discussion  of  the  physical  phenomena  and 
regimes  which  are  involved  in  nuclear  crater  formation.  For  a  more 
detailed  treatment  see  the  reviews  of  Brode*  and  Knowles  and  Brode^.  The 
bulk  of  the  energy  release  from  a  modern  nuclear  weapon  is  in  the  form  of  a 
short  (<  100  ns)  burst  of  soft  (0.5  -  2  KeV)  x-rays.  The  remainder  of  the 
energy  is  in  escaping  neutrons,  higher  energy  x-rays,  and  the  internal  and 
kinetic  energies  of  the  weapon  debris,  each  of  which  may  carry  several 
percent  of  the  total  yield.  Detailed  code  calculations  are  required  to 
determine  quantitatively  the  distribution  of  energy  among  these  forms,  but 
such  is  not  necessary  for  our  purposes. 

The  events  which  follow  a  nuclear  surface  burst  may  be  divided 
into  a  series  of  eras,  in  each  of  which  different  physical  processes  are 
dominant.  In  the  first  era  x-rays  from  the  weapon  stream  into  the  ground 
and  air,  heating  them  to  temperatures  between  10^  k  and  lO'’  K,  lowering 
("burning  out")  their  opacity,  forming  a  hemispherical  fireball  in  the  air, 
and  raising  the  pressure  in  the  surface  layers  of  the  ground  to  ~  1000 
Mbar.  After  the  weapon  has  radiated  most  of  its  internal  energy,  in 
~  100  ns,  and  lost  most  of  the  remainder  in  adiabatic  expansion,  it  ceases 
to  be  a  source  of  photons  although  the  debris  kinetic  energy  may  be 
significant.  The  heated  regions  of  air  and  ground  continue  to  expand  as 
radiation  diffuses  away  from  the  source. 
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On  times  >  10”^s  hydrodynamic  motion  becomes  important,  as  the 
heated  ground  surface  blows  off  in:o  the  air,  sending  a  shock  downward. 
Hydrodynamics  replaces  radiation  as  the  mechanism  of  energy  transport, 
first  in  the  opaque  ground,  later  in  the  transparent  air.  As  the  energy  is 
distributed  over  more  and  more  matter,  temperatures  and  pressures  drop  and 
opacities  rise.  The  ground  matter,  formerly  fully  ionized,  begins  to 
recombine.  Once  post-shock  temperatures  drop  below  10^  K  and  pressures 
below  1  Mbar,  the  ground  is  no  longer  ionized  but  only  vaporized.  At 
somewhat  lower  temperatures  and  pressures  there  is  no  vaporization.  If  the 
ground  is  wet  then  vaporization  remains  important  for  a  while  longer;  the 
nuclear  explosion  has  produced  a  steam  explosion.  But  eventually  even 
water  vaporization  stops  and  one  enters  the  regime  of  nearly  incompressible 
hydrodynamic  flow. 

In  the  era  of  hydrodynamic  flow  stresses  still  far  exceed  the 
yield  strength  of  rock  but  are  not  sufficient  to  change  significantly  its 
volume.  As  the  shock  continues  to  spherically  diverge  and  weaken,  the 
amount  of  energy  dissipated  (its  reversibility)  rapidly  drops.  Eventually 
stresses  drop  to  a  few  kilobars.  In  this  regime,  crucial  for  determining 
the  size  of  the  final  crater,  the  detailed  mechanical  properties  of  the 
ground  are  important.  Its  strength  counts,  as  does  its  behavior  once  it 
has  yielded.  It  fractures  but  is  not  completely  pulverized;  the  fractured 
material  under  pressure  has  some  strength,  like  beach  sand.  Material 
strength  and  (for  large  craters)  gravity  bring  the  cratering  process  to  a 
halt  on  a  time  scale  of  a  few  seconds.  In  some  materials  there  may  also  be 
late  time  slumping  and  other  adjustments  of  the  crater  shape.  These  too 
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involve  the  interplay  between  material  strength  and  gravity  and,  when 
present,  extend  the  cratering  process  to  minutes  or  longer. 

In  HE  (chemical)  explosions  temperatures  are  much  too  low  for 
radiative  exchange  to  be  important.  The  cratering  process  begins  with 
pressures  which  are  characteristically  100  Kbar  and  temperatures  of  the 
order  of  a  few  times  10-*  K.  HE  cratering  begins  in  the  regime  of 
incompressible  flow  and  explores  only  the  last  of  the  eras  involved  in 
nuclear  cratering.  It  is  nevertheless  of  great  interest  since  it  is 
precisely  this  last  era  which  is  so  sensitive  to  material  equation  of  state 
issues  and  most  difficult  to  treat  theoretically.  Furthermore,  once  the 
nuclear  cratering  process  enters  the  incompressible  flow  regime  it  has  lost 
all  memory  of  what  went  before.  Therefore  one  can  in  principle  design  HE 
experiments  which  reproduce  all  of  the  late  time  processes  involved  in 
nuclear  cratering. 


III  THEORETICAL  AND  PHENOMENOLOGICAL  SCALING  LAWS 
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3.1  Dimensional  Analysis  and  Centrifuge  Scaling  Laws 

We  begin  by  reviewing  the  application  of  dimensional  analysis  to 
the  problem  of  cratering.  This  has  been  discussed  in  various  degrees  of 
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completeness  by  many  authors  over  the  years.  Most  recently  Schmidt  and 

Holsapple^  have  treated  the  subject  with  considerable  care  and 
thoroughness. 

First  consider  high  explosive  (HE)  bursts  which  involve  the 
mechanical  coupling  of  gaseous  explosion  products  with  the  ground  and 
air.  We  assume  throughout  that  the  ground  is  homogenous  and  unlayered. 
Although  there  has  been  controversy  in  the  past,  it  should  be  obvious  that 
the  particular  set  of  parameters  chosen  to  describe  the  HE  source,  the 
ground,  and  the  auxiliary  variables  have  no  effect  on  the  conclusions  one 
can  draw  from  dimensional  analysis  provided  one  always  deals  with  a 
complete  set  of  variables.  For  the  problem  at  hand  we  have  three 
fundamental  dimensions:  length  (L),  mass  (M)  and  time  (T). 

For  definiteness  we  choose  variables  to  describe  the  HE  source  as 

follows : 


Source  variables 

Symbol 

Dimensions 

Total  energy  release  (yield) 

E 

[ML2/T2J 

Source  specific  energy 

Qe 

[l2/t2] 

Source  mass  density 

6 

[M/L3] 
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The  quantities  Qfc  and  6  are  intensive  variables  and  constant  for  a 
given  type  of  HE;  E  is  an  extensive  variable. 

Since  one  can  build  from  these  quantities  the  fundamental  units 
of  L  ,  M  ,  and  T  ,  it  follows  that  one  can  express  all  other  variables 
which  characterize  the  source  as  dimensionless  variables  u  ,  it  .  .  . 
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multiplied  by  appropriate  powers  of  E  ,  Qe  and  6  .  Some  of  these  other 

source  variables  can  be  expressed  trivially  in  terms  of  E  ,  Qe  and  6 

[e.g.,  source  mass  M  =  E/Qe  and  source  volume  V  =  E/ ( Qe  6)]  and  hence 

can  be  ignored  completely.  All  others  are  to  be  regarded  as  additional 

parameters  which  are  required  to  completely  specify  the  source  (e.g., 

t^X i  *  Y  »  the  specific  heat  ratio  for  the  HE  detonation  products; 

—  1  /? 

tt  =  U(Q  )  ,  where  U  is  the  characteristic  (Chapman-Jouget) 

At  G 

velocity  of  the  detonation  products;  it  =  source  aspect  ratio  (for  non- 

A  j 

spherical  sources);  etc.).  Not  all  of  these  latter  variables  are 

independent.  For  example,  if  an  ideal  gas  equation  of  state  is  assumed  for 
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the  detonation  products,  i^  =  (i^j  -  1)  .  However,  the  equation  of 

state  of  the  explosive  does  not  need  to  be  known  unless  one  seeks  to 
compare  cratering  experiments  done  with  different  types  of  HE.  We  will 
assume  nothing  about  the  equation  of  state  of  the  HE  products. 

For  nuclear  weapons  which  release  energy  largely  in  the  form  of 
radiation  a  more  natural  set  of  source  variables  might  be  total  yield  E  , 
volume  energy  density  u  -  aT^  ,  and  mass  M  along  with  appropriate 
dimensionless  variables.  Since  one  can  form  fundamental  variables  of  L  , 

M  ,  and  T  from  these  source  variables,  dimensional  analysis  for  nuclear 
cratering  proceeds  in  a  manner  identical  to  that  for  HE  cratering. 


Our  choice  of  soil  variables  is  the  following: 


Soil  variables 


Mass  density 
Yield  strength 
Viscosity  (kinematic) 


Symbol 


Dimensions 


[M/LJ] 

[M/LT2] 

[L2/T] 


All  other  variables  which  characterize  the  soil,  for  arbitrarily  complex 
equations  of  state  including,  for  example,  elastic,  plastic,  hysteretic, 
and  fracture  behavior,  can  be  expressed  as  dimensionless  variables 


"SI  ’  1IS2  * 


.  times  powers  of  the  above  three  soil  variables. 


It  should  be  emphasized  that  the  term  "viscosity  "  is  being  used 

here  in  a  generic  sense.  It  could  be  an  actual  viscosity  parameter 

describing  soil  behavior,  or  a  parameter  constructed  from  Y  ,  p  and  some 

characteristic  length  or  time  scale  of  the  medium.  For  example,  in  an 

inhomogenous  medium  made  up  of  grains  of  size  £  one  can  construct  an 

1/2 

effective  viscosity  parameter  v  ^  “  £(Y/p)  •  Similarly  heat 

conduction  introduces  a  viscosity  v  *=  k/cp  where  k  is  the  thermal 

eff 

conductivity  and  c  the  specific  heat  of  the  ground  material.  Of  course, 
once  one  viscosity  parameter  is  introduced,  all  others  can  be  accounted  for 
by  dimensionless  variables  Hg^  . 


To  illustrate  the  above  consider  the  popular  Coulomb-Mohr-Drucker- 
Prager  model  which  is  characterized  by  the  following  quantities:  elastic 
constants  K  (bulk  modulus),  o  (Poisson's  ratio),  and  constants  C 
(cohesion)  and  <f>  (angle  of  internal  friction)  which  specify  the  shear 
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value  Y  =  C  +  P  tan  $  at  which  elastic  behavior  fails  and  plastic 
behavior  commences  when  the  material  is  under  hydrostatic  pressure  P  . 

Also  required  are  rules  for  plastic  behavior  and  fracture.  The  plastic 
regime  in  this  model,  although  complicated,  is  described  by  an  equation  in 
which  strain  rate  is  proportional  to  stress  rate  (times  a  non-linear 
function  of  the  stresses).  Such  a  relation  describes  a  situation  in  which 
there  is  no  viscosity-like  p*  -ameter  for  the  soil.  Specification  of 
fracture  in  the  model  requires  an  additional  parameter  Y^  (and  perhaps 
one  or  more  dimensionless  parameters)  as  well.  For  brittle  materials  like 
rock  the  plastic  regime  is  essentially  absent;  one  can  regard  the 
parameters  C  and  $  as  specifying  directly  the  conditions  for 
fracture.  For  this  latter  case  one  could  select  as  a  complete  set  of  soil 
variables:  Y  -  C  ,  ^S1  “  K/C  >  *s2  =  0  ,  -  $  ,  and  v  =  0  .  In 

what  follows  we  will  not  restrict  our  discussion  by  assuming  any  particular 
model  of  soil  behavior. 

In  addition  to  the  source  and  soil  variables  we  must  consider  the 
following  additional  quantities: 


Auxiliary  variables 

Symbol 

Dimensit 

Gravity 

E 

[L/T2] 

Height  (depth)  of  burst 

h 

IL] 

Air  pressure 

Pa 

[M/LT2] 

Air  density 

pa 

[M/L3] 

Air  viscosity 

Ua 

[M/LT] 

We  can  form  two  independent,  dimensionless  ratios  from  this  list 


*A1  =  Pa8h/Pa 


(3.1a) 


i 


Z 


and 


a 


which  we  shall  use  in  place  of  p  and  p  respectively.  (The  special 

case  h  -*•  0  will  be  discussed  below.)  In  general  one  has  further 

dimensionless  variables  ,  n..  ,  .  .  .  which  characterize  the  air 

A3  A4 

medium. 


Our  final  task  is  to  form  dimensionless  ratios  of  soil  to  source 
variables  and  auxiliary  to  source  variables.  Since  we  have  three 
fundamental  dimensions  (  M  ,  L  ,  and  T  )  there  are  in  general  three 
independent  ratios  of  each  type.  With  no  loss  of  generality  we  choose  for 
the  soil-source  ratios: 


ft 


e  Y 

"S/Xl  '  6Q 

e 


"s/xz 


£ 

6 


ft 


> 


"s/xs 


v6 
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EL/3  1/6 
e 


and  for  the  auxiliary-source  ratios: 


(3.2) 

(3.3) 

(3.4) 


1 1 1—5 


"a/xi 


ee1/3 

Q4/351/J 

e 


(3.5) 


"A/X2  ‘  ^173 


(3.6) 


p 

"A/X3  =  13^7  *  °’7) 

We  are  now  in  a  position  to  consider  any  dependent  variable  of 
interest:  crater  volume,  crater  depth,  peak  shock  at  range  r  ,  etc.  If 

we  choose  to  consider  apparent  crater  volume,  dimensional  analysis  assures 
us  that 


v  *=  L  F  (irs/x  ,  nA/x  »  *x  ,  *s  »  nA  )  (3.8) 

where  F  is  some  (unknown)  dimensionless  function  of  its  arguments  and 
L  is  any  length  scale  in  the  problem.  (In  the  above  Ug^x  stands 
collectively  for  i?s/xl  ,  *s/x2  ,  "s/x3  ;  *x  for  irxl  ,  *x2  .  •  •  •  J 
etc.)  To  different  choices  of  L  correspond  different  functions  F  ;  all 
physical  conclusions  remain  the  same,  however,  for  any  choice  of  L  .  For 
definiteness  we  choose  L  to  be  the  linear  dimension  of  the  HE  source 
L  =  (E/6Qe)^3  and  rewrite  Eq .  (3.8)  in  the  form 


F  f  7t 


S/X  ’  A/X 


(3.9) 


where 


1 
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Equation  (3.9)  expresses  the  dependence  of  crater  volume  on  source,  soil 
and  auxiliary  variables  for  the  general  case.  Analogous  expressions  can  be 
written  down  for  all  other  dependent  observables. 


Now  let  us  consider  varying  E  for  a  fixed  HE  type,  a  fixed  soil 

type  and  a  fixed  atmosphere.  We  ask  if  it  is  possible  to  vary  other 

parameters  in  the  problem  in  a  manner  which  keeps  all  the  arguments  of  F 

fixed.  This  is  called  a  similarity  transformation.  By  hypothesis  the 

trx  ,  Hg  ,  and  ,  .  .  .  variables  are  unchanged  and  thus  can  be 

ignored.  (We  discuss  it.,  and  momentarily.)  By  inspection,  we  see 

A1  A2 

that  and  itg/x2  are  also  fixed  but  that  Tg/X3  which  involves  the 

soil  viscosity  parameter  v  is  not  fixed.  We  see  too  that  1IA/X3  is 

fixed,  1IA/X2  can  ke  kept  fixed  provided  we  scale  h  according  to 
—1/3 

hE  =  const.  ,  and  that  n  ^  WH1  remain  constant  if  we  scale  g 
according  to  gE*^  =  const.  Both  of  these  requirements  are  compatible 
with  the  required  constancy  of  n  since  the  latter  is  proportional  to 
the  product  gh  .  However,  the  dimensionless  air  viscosity  parameter 
is  not  invariant  under  these  scalings. 


We  conclude  that  similar  experiments  are  not  possible  even  if 
variations  of  g  and  h  are  allowed  whenever  the  soil  and/or  air  has  any 
significant  rate  dependent  behavior.  When  soil  and  air  viscous  behaviors 
are  absent,  however,  similarity  is  possible;  for  a  fixed  type  of  HE,  soil, 
and  atmosphere,  one  has  for  this  special  case 
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where  for  compactness  we  have  omitted  the  it  and  ir  and  any  additional 

A  D 

it  variables  in  the  argument  list  of  F  ,  since  they  are  all  fixed. 

A 


Although  Eq.  (3.10)  is  general,  it  is  inconvenient  for  discussing 
the  limit  h  +  0  since  h  appears  twice  in  the  argument  list.  It  is 
equally  general  to  retain  only  one  of  the  h  dependent  variables  and 
replace  the  other  by  the  ratio  of  the  two  h  dependent  variables 
multiplied  as  desired  by  powers  of  one  or  more  of  the  other  variables  in 
the  argument  list  of  Eq.  (3.10).  Hence  an  equivalent  to  Eq.  (3.10)  is 


3.1.1  Zero  Height  of  Burst 

If  we  continue  to  discuss  the  case  where  viscosities  play  no 
important  role  and  specialize  to  zero  height  of  burst  then  Eq.  (3.10') 
implies 


(3.11) 
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t 


or  equivalently 


\  =  G(1,S/X1  ’  "s/XZ  ’  ^A/Xl  ’  ^7X3  ’  1,/) 


(3.11') 


where  G  stands  for  F'  with  its  fourth  argument  set  to  zero  and 
Pa  ^*1  .  1'a/X3 

it'  =  -p—  =  - - -  has  been  introduced  for  convenience.  Certain 

6  nA/X2  •  "a/XI 

features  of  Eq .  (3.11)  can  be  tested  by  doing  experiments  in  accelerated 
reference  frames,  centrifuges  being  the  most  convenient.  Note  that  for 
fixed  soil,  air  and  HE  type  the  only  independent  variable  is  the 
combination  gE^^  .  Thus  either  g  or  E  or  both  may  be  varied  to 
determine  the  functional  form  of  G  . 


To  our  knowledge  the  only  relevant  centrifuge  experiments  which 
have  been  carried  out  to  date  are  those  of  Schmidt  and  Holsapple  at 
Boeing. These  complement  the  work  done  at  normal  gravity  by 
Piekutowski^  at  the  University  of  Dayton  Research  Institute.  The  early 
Boeing  work  was  done  in  dry  Ottawa  sand,  a  material  which  has  little  or  no 
cohesion  (Y  «>  0)  and  is  expected  to  have  little  or  no  rate  dependent 
effects  (v  =  0)  .  Some  results  are  shown  in  Fig.  3.1.  The  points  labeled 

"Similarity  Test"  correspond  to  two  shots  at  different  g  values  with 

1/3 

yields  differing  by  a  factor  of  3  but  having  a  common  value  of  gE  ' 

The  expected  similarity  is  confirmed  and  the  reproducibility  of  the  results 
demonstrated  as  well.  (Actually  as  indicated  on  the  figure  the  effective 
values  of  6  and  for  the  two  cases  differ  slightly  because  of  the 
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presence  of  the  initiator  which  constitutes  a  nontrivial  fraction  of  the 

explosive  for  these  rather  small  shots.  In  recognition  of  this,  the 

authors  of  ref.  5  have  cho1  to  adjust  by  selecting  g  to  keep 

n^/Xi  =  (g/Q^  j  E/6Qe  constant  rather  than  simply  •  Also  shown 

in  Fig.  3. 1  is  a  third  shot  at  a  nearby  value  of  n. which  gives  a 

A/Xi 

nearby  value  for  scaled  crater  volume  thus  indicating  similarity  over  an 

8.3:1  range  of  E  .  While  this  is  not  an  extensive  range  of  E  it  does 

give  strong  experimental  support  to  the  notion  of  centrifuge  scaling. 

Taken  seriously  these  shots,  done  with  a  few  grams  of  explosive  at 

g  =  (300  -  460)  g  ,  are  equivalent  to  TNT  yields  at  normal  gravity 
o 

(E  «  g^)  of  E  =  6  x  1018  ergs  =  0.14  KT  (gQ  =  terrestrial  value  of 

gravity). 


It  is  interesting,  too,  that  dry  sand,  although  not  expected  to 
display  actual  viscous  behavior,  does  have  a  characteristic  length  scale 
set  by  grain  size  which  as  described  above  can  act  as  an  effective 
viscosity.  The  authors  of  ref.  5  report  that  crater  sizes  at  the  g 
values  used  in  their  work  seem  to  be  insensitive  to  grain  size. 

There  is  also  another  hidden  "size"  parameter  in  the  centrifuge 
work.  Since  the  sample  size  is  finite  in  depth  and  radius  it  is  possible 
that  reflections  from  the  walls  or  bottom  of  the  sample  container  influence 
the  final  crater  size  or  shape.  It  is  easy  to  test  for  the  presence  of 
such  effects  by  varying  sample  sizes,  by  putting  absorbing  material  on  the 
sample  boundaries,  etc.  One  can  also  test  experimentally  for  other 
unwanted  influences  such  as  Coriolis  effects  and  surface  winds.  Schmidt 
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and  Holsapple  report  that  there  is  no  sign  that  these  latter  influences  are 
present  in  their  work  thus  far,  particularly  for  dry  granular  material. 


Dimensional  analysis  says  nothing  about  the  functional  dependence 

of  G  upon  its  arguments,  in  particular  upon  *  By  making 

additional  assumptions,  some  simplified  cases  can  be  examined.  If  we 

consider  the  limit  of  a  cohesionless  material  Y/ ( 6Q  )  +0  ,  negligible 

e 

atmospheric  effects  P  /<5Q  -*•  0  ,  p  / 6  -*•  0  ,  and  assume  that  the  crater 

a  e  a  - 

volume  V  is  independent  of  the  source  specific  energy  density  ,  it 
follows  that  G  must  be  of  the  form 

1/3  v-a 

G(0,  -,  u  ,0,0)  «  (u  )  8  f(^)  =  ( — — - — \  8  f(-) 

6  A/ XI  J  1  A/X1J  V  I  4/3  1/3  |  V 


(Q  6 

e 


(3.12) 


where  a  =  3/4  since  Qe  appears  linearly  on  the  left  hand  side  of 
S 

Eq.  (3.11).  Here  f  is  an  unknown  function  of  its  argument  but  is 
constant  for  fixed  HE  and  soil  combinations.  Hence  for  this  special  case, 


„  rE  x3/4 

V  =  const  I — J 

g<5 


(3.13) 


a  result  often  referred  to  as  gravity  scaling. 

One  can  ask  if  Eq.  (3.12),  or  equivalently  (3.13),  is  satisfied  by 
HE  bursts.  The  answer  for  Ottawa  sand  is  negative  as  shown  in  Fig.  3.1. 
Remarkably  the  functional  dependence  of  G  for  this  material  over  a  wide 
range  is,  however,  well  fitted  by  the  simple  form 
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const  (*A/X1)a  , 


(3.14) 


it 


V 


where  the  experimental  value  of  the  exponent  is  a  «  0.472  ±  0.010  , 

clearly  different  from  a  =  3/4  .  This  difference  means  that  the  crater 

g 

volume  for  this  material  depends  not  only  on  the  total  energy  release  E 
but  also  on  the  specific  energy  Qg  •  This  was  also  directly  demonstrated 
by  comparison  of  lead  azide  and  PETN  shots.  An  equivalent  statement  is 
that  V  depends  on  E  and  the  source  volume  =  (E/ 6Qe )  ,  i.e.,  the 
crater  "remembers"  the  size  of  the  explosive.  This  violates  no  general 
principles  but  shows  that  the  key  assumption  behind  the  gravity  scaling 
argument  is  incorrect  for  values  of  Qe  in  the  range 
(1-5  x  lO*®  ergS/gm)  used  in  the  experiments.  The  dependence  or 
independence  of  V  on  Qg  in  other  ranges  is  an  open  question  but  clearly 
an  important  one,  when  one  contemplates  extrapolation  to  nuclear  energy 
densities.  Hypervelocity  impact  cratering  is  therefore  of  great  interest 
since  it  extends  Qg  to  values  well  above  the  range  accessible  with 
chemical  explosives.  In  detailed  experiments,  Schmidt  and  Holsapple  have 
observed  that  the  dependence  expressed  by  Eq .  (3.14)  holds  also  for 

impact  velocities  from  0.5  km/s  to  7  km/s  which  corresponds  to  a  200:1 
ratio  for  Qg  and  a  maximum  Qg  of  2  x  10**  ergs/gm. 


A  second  special  case  is  the  limit  in  which  gravity  effects  are 
negligible  compared  to  material  strength  effects.  In  this  case  the 
dimensionless  variable  it^^j  drops  out  of  the  problem  and  Eq.  (3.11) 
reduces  to 
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(3.15) 


This  is  the  familiar  "cube  root"  scaling  law.  Although  this  law  is  often 
assumed  in  crater  scaling  phenomenology  it  has  no  fundamental  basis.  In 
contrast,  for  air  bursts  it  is  physically  plausible  that  gravity  plays  an 

Q 

insignificant  role;  indeed,  cube  root  scaling  is  observed  as  predicted.0 
Note  that  even  when  it  applies,  Eq.  (3.15)  says  nothing  about  the 
dependence  of  crater  volume  on  material  properties.  To  determine  this  one 
must  do  actual  or  numerical  experiments. 

Centrifuge  experiments^  in  Kirtland  Air  Force  Base  (KAFB)  alluvium 
exhibit  an  approximate  cube  root  regime  as  shown  in  Fig.  3.2  and  also  a 
nontrivial  dependence  on  explosive  type.  Fig.  3.2  also  shows  that  at  high 
g  values  crater  volumes  in  alluvium  and  clay  no  longer  grow  linearly  with 
energy  release.  Holsapple  and  Schmidt^  have  been  able  to  find  a  simple 
parametrizatlon  of  all  of  their  data  in  terms  of  two  ( Coulomb-Mohr) 
material  parameters:  cohesion  and  angle  of  internal  friction 


V 


(0.174) 


{tan  $  +  0. 1 } 


-0.472 


(3.16) 


Curves  corresponding  to  this  formula  are  shown  on  Fig.  3.2  and  represent 
well  the  general  features  of  the  data.  Centrifuge  results  for  wet  sands'* 
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are  also  describable  by  Eq .  (3.16)  using  an  effective  value  of  q,  less 
than  that  for  dry  sand  and  a  non-zero  effective  cohesion  parameter  C  . 

It  is  interesting,  too,  that  the  above  formula  based  only  on 
centrifuge  studies  correctly  predicts  the  apparent  crater  volumes  obtained 
in  alluvium  for  TNT  bursts  of  5000  and  256  lbs  done  in  the  field  at  the 
Nevada  Test  Site  (NTS) — see  Fig.  3.3.  (Note  that  these  data  are  only 
barely  into  the  gravity-influenced  regime.  The  correction  factor  to  cube 
root  scaling  due  to  the  gravity  term  in  Eq .  (3.16)  corresponds  to  a 
reduction  in  apparent  crater  volume  of  only  11%  and  25%  for  the  lighter  and 
heavier  shots  respectively.) 

Although  the  form  of  Eq.  (3.16)  can  be  partially  motivated  by 
energy  balance  considerations,  it  must  in  the  end  be  regarded  as  a 
phenomenological  expression  which  may  require  modification  for  materials 
more  complex  than  those  used  in  its  development.  It  does,  however,  stand 
as  an  important  reminder  that  there  can  be  various  regimes  (cohesion 
dominated,  gravity  influenced,  gravity  dominated,  etc.)  in  cratering  which 
render  naive  extrapolations  meaningless. 

The  dependence  of  crater  radius  r  and  depth  d  on  material 
properties  and  source  variables  has  also  been  explored  in  the  centrifuge 
crater  work.^0  The  patterns  seen  in  the  radius  and  depth  variables  are 
somewhat  more  complicated  than  those  described  by  the  volume  formula, 

Eq.  (3.16).  As  explosive  yield  moves  from  the  cohesion  dominated  regime 
(small  E  at  gQ  )  into  the  gravity  influenced  regime  (larger  E  at  g0  ) 
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Figure  3.3  CRATER  VOLUME  VS  YIELD  (EQUIVALENTLY  MASS!  FOR  TNT  AT  ZERO 
HEIGHT  OF  BURST  OVER  ALLUVIUM.  THE  CURVE  CORRESPONDS  TO  EQ.  (3.161  WITH 
PARAMETERS  DETERMINED  ENTIRELY  BY  CENTRIFUGE  CRATER  STUDIES.  (Source; 


the  crater  aspect  ratio  r/d  increases  suggesting  that  stability 
considerations  are  becoming  important.  The  same  trend  of  increasing  aspect 
ratio  with  increasing  crater  size  is  well  known  in  the  PPG  nuclear  crater 
results  (see  Fig.  3.4)  and  also  in  earth  impact  crater  observations.  Of 
course  in  these  latter  two  cases  geological  layering  and  other  phenomena 
may  also  be  playing  a  role.  Nevertheless,  the  trend  of  the  centrifuge 
results  are  notable  and  would  seem  to  indicate  that  phenomena  important  to 
large  scale  craters  are  being  probed  by  the  centrifuge  studies. 


3.1.2  Non-Zero  Heights  of  Burst 

All  of  the  above  centrifuge  crater  results  were  for  zero  height  of 
burst  which  by  definition  correspond  to  a  half  buried  spherical  charge 
distribution.  Only  very  limited  results  are  presently  available  for  non¬ 
zero  h  values.  Most  are  for  tangent  above  charge  emplacements,  i.e., 
h  =  +a  (charge  radius)  over  Ottawa  sand.  Preliminary  results1*  indicate 
that  for  this  cohesionless  material  crater  volume  departs  even  more 
strongly  from  cube  root  scaling  for  positive  h  than  the  fit  to  the  h  =  0 
data  described  by  Eq.  (3.14).  There  seems  to  be  an  additional  dependence 


of  u  on  source  specific  energy  beyond  that  included  in  the  dimensionless 
variables  (gravity  scaled  yield)  and  nA/x2  (scaled  height).  This 
could  indicate  that  the  scaled  atmospheric  pressure  variable  1T^/x3  ls 


playing  a  role. 


For  negative  values  of  h  Schmidt  and  Holsapple  report11  that  two 
tangent-below  (h  =  -a)  shots  of  the  same  h  type  in  Ottawa  sand  at 
different  g  and  E  values,  but  common  gE  '  ,  show  similarity  as  was 


I 
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Figure  3.4  THE  TREND  OF  CRATER  ASPECT  RATIO  R/D  WITH  CRATER  DEPTH  D  FOR  SURFACE 

AND  NEAR  SURFACE  BURSTS  FOR  NUCLEAR  AND  HE  SOURCES  AT  NORMAL  GRAVITY 
AND  HE  CRATERS  AT  ENHANCED  GRAVITY.  THE  CENTRIFUGE  CRATER  DEPTHS  HAVE 
BEEN  SCALED  BY  g  AS  REQUIRED  BY  SIMILARITY.  THE  CENTRIFUGE  CRATERING 
MEDIUM  WAS  WET  SAND.  (Source:  ref.  9) 


observed  at  h  =  0  .  These  workers  have  also  achieved^  successful 
centrifuge  simulations  in  alluvium  of  the  20  T  Stagecoach  III  (h  =  -3.77  a) 
and  20  T  Jangle  HE  2  (h  =  -1.19  a)  TNT  shots.  These  complement  the 
successful  (±  20%)  simulation^  of  the  shallow  buried  0.5  KT  nuclear  Johnie 
Boy  event  (h  »  -0.6A  a)  done  early  in  the  centrifuge  program.  The 
centrifuge  HE  charge  equivalent  for  Johnie  Boy  was  established  by  applying 
the  centrifuge  scaling  law  gE  =  const,  to  the  full  scale  HE 
equivalent  charge  calculated  in  ref.  (13). 

3. 2  Phenomenological  Scaling  Laws 

It  has  not  been  possible  to  obtain  scaling  laws  for  near-surface 
burst  cratering  by  a  direct  study  of  the  underlying  equations  of  motion  and 
equations  of  state.  This  contrasts  with  the  rather  well  understood 
situation  for  the  blast  wave  produced  by  atmospheric  bursts.  Consequently 
a  great  deal  of  work  has  been  devoted  to  the  development  of 
phenomenological  scaling  laws  based  on  field  data.  Many  HE  cratering 
experiments  have  been  done  co  further  these  efforts.  A  recent  and 
comprehensive  scaling  phenomenology  is  that  developed  by  Cooper^  and 
coworkers*-*;  it  forms  the  basis  of  the  prediction  curves  given  in  the  Air 
Force  Manual  for  Design  and  Analysis  of  Hardened  Structures^  and  in  the 
DNA  Effects  Manual.^  Alternative  phenomenologies  exist  as  well,  e.g. , 
Refs.  (3),  (A)  and  (I7a).  The  latter  approach^3  due  to  Lilley  challenges 
the  notion  that  nuclear  sources  must  be  separated  into  low  energy  density 
and  high  energy  density  categories  as  done  by  Cooper,  et  al. 

Strong  assumptions  must  be  made  concerning  the  form  of  the 
phenomenological  scaling  laws  in  order  to  make  any  progress.  There  do  not 
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exist  sufficient  data  to  fully  test  these  assumptions.  Hie  primary 
assumption  of  Cooper  is  that  the  apparent  crater  volume  V  is  given  by  the 
product  of  three  factors:  yield  E  (cube  root  scaling),  a  constant  n 
which  contains  all  the  effects  due  to  material  properties  and  geology,  and 
p  ,  a  function  which  depends  on  the  scaled  depth  of  burst  D  and  weapon 
type  i  [i  =  HE,  h’E  (non-radiating),  NE  (radiating)].  Namely 


V  =  En<f>.(D) 


(3.17) 


With  no  loss  of  generality  we  can  take  $HE  (0)  =  1  .  Then  by  virtue  of 
the  assumed  form  of  Eq.  (3.17),  the  values  of  (0)  and  <t>fjgn_ra^  (0) 

give  the  so-called  relative  efficiencies  for  cratering  by  radiating  and 
non-radiating  nuclear  bursts,  respectively,  in  any  geology.  Similarly  the 
ratio  ra/g  -  h(a)/  n(6)  gives  the  relative  cratering  efficiency  in 
geology  a  to  that  in  geology  g  for  any  type  source  assuming  the 
validity  of  Eq.  (3.17). 


Central  to  Cooper's  approach  is  the  further  assumption  that  all 

1/3 

dependent  and  independent  lengths  are  to  be  scaled  by  V  .  In 
particular  in  Eq.  (3.17) 


D  =  d/V1/3 


(3.18) 


where  d  is  the  actual  depth  of  burst.  (Note  d  here  corresponds  to  -h 
used  in  subsection  3.1  above.)  As  a  consequence  of  this  scaling  rule  the 
dependent  variable  V  appears  on  both  sides  of  Eq.  (3.17).  However,  by 
rewriting  Eq.  (3.17)  in  the  form 
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(3.19) 


and  solving  for 


V_ 

nE 


we  can  re-express  it  in  the  form 


V  = 


(3.20) 


This  is  a  canonical  cube  root  scaling  law  in  terms  of  an  "effective"  yield 
nE  .  Note  that  the  height  (depth)  of  burst  function  depends  on 

source  type  i  but  is  independent  of  material  properties;  all  dependence 
on  the  latter  has  been  swept  into  an  effective  yield  variable. 


Equation  (3.20),  or  equivalently  (3.17),  implies  strong  statements 
about  the  crater  volume.  In  particular  the  cube  root  dependence  allows  no 
dependence  on  gravity  for  even  the  largest  craters.  Evidence  for  such  a 
cube  root  scaling  law  for  HE  surface  bursts  seems  to  be  based  mostly  on 
field  tests  done  with  256  and  1000  lb  TNT  charges  in  cohesive  soils 
(alluvium  tuff,  playa,  clays,  wet  sands,  etc.).  Reproducibility  in  a  given 
geology  shows  scatter  of  a  factor  of  2  to  3  due  presumably  to  site 
inhomogeneities  and  variations.10  These  variations  can  mask  small  but  real 
effects  which  might  be  important  for  extrapolations  to  yields  in  the 
kiloton  and  higher  ranges.  Indeed,  the  centrifuge  cratering  experiments 
mentioned  earlier,  when  scaled  to  normal  gravity  predict  that  departures 
from  cube  root  scaling  will  show  up  in  precisely  this  yield  range,  e.g., 
they  indicate  a  factor  of  2  effect  in  alluvium  at  E  »  1  KT  and  zero 
height  of  burst.  (See  Fig.  3.4). 
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Cooper's  fit  to  the  A  (D)  for  -0.8  <  D  <  0.8  is  shown  in  Fig. 
3.5  along  with  p  values  (cratering  efficiency)  for  generic  geologies. 

Note  there  is  considerable  spread  in  the  data  about  the  best  fit  line. 

Next  consider  near  surface  nuclear  cratering.  Here  the  database 

♦ 

(See  Table  I)  is  extremely  limited  and  field  conditions  poorly 
documented.  In  the  unclassified  literature,  there  exist  data  for  two  low 
yield,  low  energy  density  shots  done  over  alluvium  at  the  NTS  and  ten  high 
yield  shots  over  coral  at  the  PPG.  Of  these  latter,  two  (KOA  and  SEMINOLE) 
were  done  in  water  tanks  and  hence  are  usually  considered  to  be  equivalent 
to  low  energy  density  nuclear  shots.  The  remaining  eight  PPG  bursts  are 
thought  to  resemble  modern  high  energy  density  weapons  for  which  most  of 
the  energy  release  occurs  as  X  radiation.  However,  as  has  been  emphasized 
by  many  authors,  it  is  not  clear  that  any  of  these  weapons  coupled  a 
significant  amount  of  energy  to  the  ground  directly  by  radiation  since  the 
weapons  were  surrounded  by  experimental  equipment  and  sheds  and  in  some 
cases  were  detonated  on  barges. 

Because  of  the  limitations  of  the  nuclear  database  one  has  little 
chance  of  testing  the  cube  root  scaling  assumption  and  the  rule  for 

scaling  the  height  of  burst.  All  one  can  do  is  make  these  assumptions  and 
then  use  the  data  to  find  the  height  of  burst  curve  The  results 

are  shown  in  Fig.  3.6.  In  the  Cooper  model  two  separate  curves  are 

•  *• 

necessary,  one  for  low  energy  density  sources,  ij>^n  a  ,  and  one  for  high 
energy  density,  •  for  the  reasons  stated  above,  KOA  and  SEMINOLE 

have  been  combined  with  the  NTS  alluvium  shots  using  the  results  from  small 
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Figure  3.5  THE  HEIGHT  OF  BURST  CURVE  FOR  HE  (ALL  GEOLOGIES!  NORMALIZED 
ACCORDING  TO  f»HE(0>  -  1.  ALSO  SHOWN  ARE  VALUES  OF  THE  "HE 
CRATERING  EFFICIENCY"  1  FOR  GENERIC  SOURCES.  (Source:  ref.  14) 
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Figure  3.6  THE  HEIGHT  OF  BURST  CURVES  FOR  RADIATING  (HIGH  ENERGY  DENSITY) 
AND  NON-RADIATING  (LOW  ENERGY  DENSITY)  NUCLEAR  SOURCES.  THE 
VALUES  OF  THESE  CURVES  AT  D«0  GIVE  THE  RELATIVE  CRATERING  EFFI¬ 
CIENCIES  FOR  THESE  SOURCES  VS  HE  SOURCES.  THE  PRODUCT  OF  $NE(0) 
AND  THE  TJ  VALUE  FROM  FIG.  3.6  GIVES  THE  NE  CRATERING  EFFICIENCY 
ACCORDING  TO  THE  PHENOMENOLOGY  BEING  DESCRIBED.  (Adapted  from  ref.  14) 


III-26 


> 


♦ 


HE  cratering  shots  ri  ,/n  , ,  .  =  4  .  Since  the  two  NTS  and  these  two 

coral  alluvium 

PPG  craters  have  very  different  aspect  ratios,  it  is  not  clear  whether  this 
manner  of  comparison  has  any  meaning.  The  manner  in  which  the  and 

^non-rad  curves  merge  at  large  positive  and  negative  values  of  D  is  a 
matter  of  speculation  and  is  at  present  unconstrained  by  experiment. 


Note  that  when  treated  in  the  above  manner  the  dependence  of 
crater  volume,  and  presumably  other  observables,  is  a  very  steep  function 
of  height/depth  of  burst  in  the  neighborhood  of  D  «  0  .  This  is  very 
unfortunate  since  one  is  interested  in  precisely  this  near-surface  burst 
regime.  The  extreme  weakness  of  any  phenomenology  based  on  the  PPG  craters 
is  underlined  by  the  progressive  flattening  of  the  craters  (increase  in 
aspect  ratio)  with  yield  value,  a  feature  not  observed  in  either  the  HE  or 
NE  craters  at  the  NTS  nor  the  1000  lb  HE  craters  at  the  PPG — see  Fig.  3.7. 


X 
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The  very  thorough  and  extensive  restudy  of  the  PPG  craters  by 
Ristvet  et  al. ^  has  concluded  that  for  at  least  MIKE,  KOA  and  OAK  much  of 
the  apparent  crater  volume  resulted  from  compaction  of  the  underlying  coral 
medium  along  with  shape  changes  due  to  liquifaction  and  late  time  gravity 
slumping.  Extrapolation  of  PPG  data  to  other  geologies  by  simple  methods 
may  therefore  have  no  meaning  whatsoever. 


The  phenomenology  of  Cooper  and  coworkers  goes  beyond  prediction 

for  apparent  crater  volume  to  cover  such  observables  as  crater  radius, 

1/3 

crater  depth,  peak  shock  at  range  r  ,  etc.  In  all  of  these  the  V 
length  scaling  rule  is  adopted.  Since  in  all  of  these  cases  one  is  scaling 
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Figure  3.7  AVERAGED  CROSS  SECTIONS  THROUGH  HE  AND  HE  CRATERS  FORMED 
IN  PPG  CORAL.  NOTE  THE  PROGRESSIVE  FLATTENING  (INCREASE  IN  ASPECT  RATIO) 
AS  CRATER  SIZE  INCREASES.  (Source:  ref.  19) 
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a  dependent  variable  by  another  dependent  variable,  (crater  volume),  there 
is  some  reason  to  hope  that  even  if  the  assumptions  behind  Eq.  (3.17)  for 
V  are  not  fully  correct  for  nuclear  and  HE  bursts  in  the  KT-MT  range,  the 
predictions  for  other  effects  may  nevertheless  have  some  validity. 
Centrifuge  studies  could  shed  light  on  scaling.  Insufficient  data 

■» 

are  presently  available. 


3. 3  Comparison  Between  Theory  and  Experiment 

There  exists  no  case  for  which  a  computer  code  calculation  has 
been  carried  all  the  way  to  cratering  in  a  geology  and  configuration 
appropriate  to  a  PPG  radiating  nuclear  burst.  Consequently  all  comparisons 
between  code  results  and  the  PPG  database  make  use  of  the  scaling 
phemenology  discussed  in  the  last  subsection,  or  some  other  phenomenology, 
and  are  therefore  correspondingly  uncertain.  For  non-radiating  nuclear 
bursts  two  code  calculations  exist  which  go  all  the  way  to  late  state 
crater  formation:  (i)  a  calculation^  for  KOA,  a  1.3  MT  burst  exploded  in 
a  water  tank  at  the  PPG  and  (ii)  a  calculation^ 1  for  JOHNIE  BOY,  a  0.5  KT 
explosion  buried  in  NTS  alluvium  at  a  depth  of  0.585  m. 
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For  HE  bursts  the  situation  is  somewhat  better  although  the  number 
of  cases  in  which  one  can  make  actual  comparisons  of  code  results  and  field 
data  is  still  quite  limited.  Accurate  characterization  of  field  geology 
and  material  properties  has  proven  to  be  quite  complex  at  even  the  best  of 
sites.  The  best  documented  HE  crater  calculations  seem  to  be:  (i)  MIXED 
COMPANY  -  III  (0.5  KT),  a  tangent-above  burst  on  siltstone/  sandstone^ ; 
(ii)  MIDDLE  GUST  -  I  (0.02  KT),  half  buried  in  layered  clay/shale^;  and 


ie 
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(iii)  MIDDLE  GUST  -  III,  and  -  IV  (0.1  KT),  tangent-above  over  layered 
23 

clay/shale . 

Comparisons  between  code  predictions  and  experiment  for  non¬ 
radiating  nuclear  and  HE  are  shown  in  Tables  II  and  III  respectively.  One 
notes  that  the  situation  is  in  reasonably  good  shape  for  the  HE  case. 

There  is  even  good  agreement  in  some  cases  between  code  predictions  for  the 
ground  shock  patterns  and  observations  provided  material  parameters 
appropriate  to  the  in  situ  case  are  used.  This  agreement  has  only  been 
achieved  after  considerable  post  shot  refinement  of  the  computer  models, 
however. 


For  the  low  energy  density  nuclear  case  the  situation  is  mixed. 

For  the  one  small  yield  NTS  burst  the  agreement  in  all  crater  parameters  is 
quite  satisfactory  but  for  the  large  yield  PPG  burst  the  predicted  crater 
aspect  ratio  is  in  disagreement  with  the  observed  value  even  though  the 
volume  is  comparable.  Details  of  the  KOA  code  calculation^®  are  unknown  to 
us  at  the  present  time.  In  particluar  we  do  not  know  if  the  coral 
compaction  which  contributes  to  much  of  KOA  crater  volume,  as  indicated  by 
the  resurvey*'*,  was  modeled  in  the  code  calculations.  Nor  do  we  know  if 
liquifaction  and  gravity  slumping  (mechanisms  often  mentioned  as 
responsible  for  the  high  aspect  ratios  of  the  PPG  craters)  were  permitted 
by  the  code. 

As  mentioned  already,  no  direct  tests  of  the  code  calculations  for 
crater  formation  by  radiating  nuclear  weapons  are  available.  In  Table  IV 
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we  show  instead  code  predictions  vs  empirical  predictions.  The  results  in 
the  table  correspond  to  a  5.0  MT  burst  detonated  at  a  height  of  1.0  m  over 
two  layered  geologies  which  are  intended  to  be  representations  of  Minuteman 


Missile  sites.  All  three  code  calculations  share  a  common  early  time 
history  since  they  begin  with  the  output  from  a  combined  radiation- 
hydrodynamic  code  calculation  for  a  model  weapon  (Source  3).  The  ref.  25 
and  26  calculations  in  Table  IV  correspond  to  identical  geologies  with  the 
difference  that  the  material  parameters  used  in  the  first  were  modified  in 
such  a  way  that  once-shocked  material  was  assigned  a  reduced  tensile 
strength. 


The  empirical  values  shown  in  Table  IV  are  "hypothetical  data" 
constructed  by  means  of  the  Cooper  phenomenology  described  above.  Namely, 
they  involve  the  assumption  of  cube  root  scaling,  the  height-of-burst  curve 
shown  in  Fig.  3.6  and  the  geology/material  adjustments  as  determined  from 
craters  formed  by  small  HE  bursts.  Two  sets  of  empirical  radius  and  depth 
values  are  shown  in  Table  IV.  The  upper  line  corresponds  to  those 
appropriate  to  a  "bowl-shaped"  crater  and  the  lower  to  a  "disk-shaped" 
crater  as  observed  at  the  PPG.  The  upper  values  are  obtained  by 
extrapolating  HE  radius/depth  curves  and  the  lower  by  drawing  smooth  curves 
through  PPG  data.  One  notes  that  the  code  predictions  are  in  poor 
agreement  with  the  "empirical  data"  and  that  the  two  sets  of  empirical  data 
differ  significantly  as  well. 

Not  surprisingly  the  crater  is  larger  in  the  weakened  material 
(ref.  26  vs  ref.  25) — but  not  much  larger.  More  interesting  and  perhaps 
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even  encouraging  is  the  larger  aspect  ratio  observed  in  the  weakened  case, 
a  result  which  underlines  the  critical  importance  of  the  strength  of 
materials  and  fracture  criteria  employed  in  the  late  stage  of  cratering. 

Although  much  has  been  said  about  the  apparent  inability  of  codes 
to  predict  PPG-like  craters,  no  definitive  conclusions  can  be  reached.  A 
detailed  review  of  the  calculations  shown  in  Table  IV  has  been  made  by 
Kreyenhagen28  an<j  a  number  of  problems  uncovered,  especially  at  the 
Eulerian-Lagrangian  boundaries  which  appear  in  the  codes.  We  have  no 
original  insights  but  find  it  hard  to  believe  that  the  complexities  of  the 
PPG  coral  and  coral  sand  geology,  especially  the  compaction  observed  in  the 
resurveys,  can  all  be  accommodated  by  a  single  multiplicative  parameter  as 
is  assumed  when  one  uses  the  phenomenology  of  subsection  3.2  to  extrapolate 
the  PPG  observations  to  the  code  geologies.  The  simple  rules  which  work 
for  relatively  low  yield  HE  shots  need  not  work  for  large  yield  bursts. 


IV  RADIATIVE  COUPLING  AND  INITIAL  IMPULSE 


Here  we  make  a  unified,  albeit  crude,  calculation  that  includes 
both  the  inner  few  meters  that  are  directly  illuminated  by  kilovolt  X-ray 
flux  from  the  disassembling  warhead  and  also  the  next  hundred  or  so  meters 
in  radius  that  are  illuminated  by  radiation  scattered  in  the  fireball.  Our 
purposes  are  (i)  to  illustrate  the  relevant  physical  processes  in  each  of 
the  two  regimes,  and  (ii)  to  estimate,  independently  of  other 
investigations,  whether  fireball  coupling  may  be  important. 

Because  thermal  penetration  distances  are  always  very  small 
compared  with  the  scale  of  subsequent  ground  motions,  there  is  a  fairly 
clean  separation  between  "hot"  (T  >  10  eV)  and  "cold"  (T  <  1  eV) 
processes.  Heating  of  the  ground  and  its  hydrodynamic  blowoff  in  the 
fireball  is  "hot."  Groundshock,  elasto-plastic  response  of  the  rock,  and 
crater-related  ground  motions  are  "cold."  The  "hot"  processes  are  much 
more  rapid  in  time  than  the  "cold"  ones.  In  first  approximation,  the 
coupling  between  the  hot  and  cold  regimes  is  therefore  entirely  determined 
by  the  impulse  per  area  (as  a  function  of  radius)  delivered  to  the  cold 
ground.  Zeldovich  and  Raizer  (ch.  XII,  sec.  13)^  discuss  this 
approximation.  In  the  next  approximation,  the  duration  of  this  impulse 
might  also  be  relevant,  but  only  when  it  is  longer  than  the  natural 
duration  of  the  groundshock. 
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4. 1  Souxce  Model 

We  imagine  that  each  radius  r  is  illuminated  for  a  time  t(r) 
by  a  blackbody  flux  of  temperature  T(r)  .  The  functions  t(r)  and 
T(r)  are  chosen  to  represent  a  nominal  E  =  1  MT  surface  burst  at  a 
height  L  =  1  m.  For  r  <  L  ,  t(r)  is  constant: 


t 


tto(T  )4L2 
m 


(r  <  L) 


(4.1) 


Here  0/4n  «  1/4  is  “he  solid  angle  subtended  by  the  zone  of  radius  L  , 
o  is  the  Stephan-Boltzman  constant,  and  Tm  is  the  effective  (not  color) 
temperature  of  the  flux.  We  adopt  Tm  =  1.5  keV  as  a  standard  value,  and 
with  the  other  values  already  indicated,  obtain 


t(r)  =  7  x  10  ®s  *  70  ns  .  (r  <  1  m)  .  (4.2) 

For  r  >  1  m,  we  match  onto  a  "standard"  fireball  given  in  Figure  3 
and  Table  1  of  Brode*.  This  becomes  our  standard  source  and  has  t(r) 
and  T(r)  as  shown  in  Fig.  4.1. 

4. 2  Ground  Opacity  and  Radiative  Transfer 

As  an  idealization  (that  is  not  too  inaccurate)  we  take  Kramer's 
opacity  law 


K  «=  K(p,T) 


^  (k)or 


T  <  T, 


B  * 


(4.3) 


Here  Koe  is  the  electron  scattering  opacity,  K  ■  0.2  cm  /g  , 
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Figure  4.1  THE  TIME  AND  TEMPERATURE  FUNCTIONS  FOR  THE  STANDARD  1MT 
NUCLEAR  SOURCE  USED  IN  THIS  SECTION. 


We  could  be  somewhat  more  precise  and  distinguish  the  radiation 
color  temperature  from  the  ground  temperature  TG  by  using  an 

effective  opacity 


K  =  K 


i  v/T  /  T  i  J 

|£_ \|  _!  1  ( I 
Ip  )  \  td  )  It  / 

Vo/’  R/  V  q! 


(4. A) 


This  would  predict  a  somewhat  greater  penetration  of  the  radiation  in  the 
direct-coupled  inner  region.  We  have  chosen  not  to  introduce  this  extra 
complication. 


If  radiation  and  ground  are  assumed  to  have  the  same  temperature, 
then  on  scales  larger  than  the  radiation  mean-f ree-path  X  the  diffusion 
of  radiation  into  the  ground  is  governed  by 


(4.5) 


where  a  is  the  radiation  constant.  It  is  convenient  to  rewrite  Eq .  (4.5) 
in  the  form 


4acT^pr 


3K(Ts,P0)CvP: 


(A. 6) 


In  this  second  form  of  the  equation  we  have  grouped  the  dimensional 
quantities  into  a  single  diffusion  coefficient 


(4.7) 


4acT3p 

k  =  - 2-2__ 

3K(Ts^o^Cvp 

with  dimensions  (length) ^/(time) .  Here  Tg  is  the  temperature  applied  to 
the  surface,  namely  T(r)  above,  and  «  6.2  x  10^erg/g  K  is  the 

specific  heat  of  the  ground  matter. 

A  remark  on  the  dependence  of  k  on  p  is  in  order  at  this 

point.  Two  powers  of  p  appear  because  we  are  measuring  depth  in 

centimeters  rather  than  grams.  The  one  additional  power  of  p  comes  from 

Eq.  (4.3).  So  we  will  have  to  ask  (see  below)  whether  there  is  appreciable 

blowoff  of  matter  and  reduction  of  density  during  the  time  that  the 

source  T  is  applied.  Close  in,  in  the  kilovolt  region,  there  will  be 
s 

insufficient  time  for  hydrodynamic  motion,  so  we  can  hold  p  fixed  at 

p  As  a  simple  estimate  (i.e.,  which  underestimates  the  coupling  in  the 
o  - 

fireball  region)  we  can  take  p  to  be  fixed  at  its  initial  value 
everywhere . 


Treating  the  density  in  this  way,  we  see  that  the  depth  of 
penetration  of  the  half-maximum  point  of  energy  density  is  roughly 
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Hence  the  energy  deposition  per  cm 


2 


is  roughly 
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and  the  Impulse  per  cmz  generated  by  the  blowoff  of  the  photon-heated  layer 
is  roughly 

3.75 


r 

1/2 

/  T  \ 

p'  « 
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(4.10) 


In  Figure  4.2  we  plot  the  function  £(r)  as  functions  of  radius 
r  ,  using  T(r)  and  t(r)  from  our  "standard"  source. 

It  is  interesting  to  note  the  appearance  of  secondary  maxima  in 
l  at  radii  of  50  m  or  so.  This  maximum  occurs  because  the  fireball 
stagnates  for  a  long  time  prior  to  shock  breakaway.  We  now  want  to  see 
whether  this  stagnation  causes  appreciable  energy  or  momentum  deposition 
from  the  fireball  to  the  ground  at  these  large  radii. 

4.3  Energy  and  Impulse  Coupled  to  the  Ground 

Figure  4.3  plots  the  energy  (solid  curves)  and  blow-off  momentum 
(dashed  curves)  implied  by  Eqs.  (4.9  and  4.10),  with  the  T(r)  and  t(r) 
of  our  standard  source.  We  plot  r^E'  and  r^p'  vs  log  r  so  that  the 
height  of  the  curve  indicates  the  marginal  integrated  contribution  from 
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Figure  4,3  ENERGY  AND  MOMENTUM  COUPLED  TO  THE  GROUND  ACCORDING  TO  THE 
FIREBALL  COUPLING  MOOEL  OF  THIS  SECTION.  CURVES  FOR  THE  INTEGRATED  CON¬ 
TRIBUTIONS  WITH  AND  WITHOUT  THE  QUESTIONABLE  INTERMEDIATE  REGION  ARE 
SHOWN. 


dr  .  Also  plotted  for  both  energy  and  momentum  is  an  integral  curve 
showing  the  total  energy  (in  kiloton  equivalent)  and  momentum  (in  cgs 
units)  deposited  inside  radius  r  .  In  each  case,  we  have  computed  the 
integral  in  two  different  ways,  either  including  or  excluding  the 
contribution  predicted  by  Eq.  (4.9)  or  Eq .  (4.10)  between  2  m  and  20  m 
radius.  This  point  requires  some  discussion: 

Consider  the  curve  of  integrated  energy  deposition.  Inside  1  m 
radius,  our  calculation  predicts  35  KT  (i.e.,  3.5%  efficiency)  deposited 
into  the  ground.  This  number  is  in  line  with  the  detailed  calculations^? , 30^ 
as  is  our  value  for  the  central  penetration  depth  of  the  radiation. 

Between  1  m  and  2  m,  the  integral  rises  to  about  9%  energy  efficiency, 
which  is  not  too  far  wrong,  considering  the  crudeness  of  our  treatment. 
Between  2  m  and  20  m,  however,  it  rises  to  more  than  50%.  This  is 
certainly  spurious,  and  occurs  for  a  simple  reason:  the  model  assumes  that 
the  fireball  illuminates  the  ground  as  a  black  body,  i.e.,  that  the 
radiation  field  within  the  fireball  is  spatially  isotropic.  But  the 
scattering  length  of  radiation  in  the  fireball  is  of  order  tens  of 
meters.  Therefore,  until  the  fireball  is  of  order  20  m,  the  free-streaming 
radiation  field  does  not  (for  the  most  part)  illuminate  the  ground 
isotropically.  Our  effective  temperature  idealization  is  fairly  good  in 
the  first  meter  or  so,  where  the  illumination  is  from  the  source  directly; 
it  ought  also  to  be  good  outside  of  20  m,  where  the  fireball  is  optically 
thick  to  scattering,  but  should  not  be  used  in  between.  Ihe  choice  of  20  m 
as  a  cutoff  radius  is  not  critical,  since  the  curves  in  Figure  4.3  show 
very  little  additional  contribution  between  15  m  and  40  m. 


Between  40  m  and  100  m,  there  is  a  significant  additional 


deposition,  however.  Integrated  energy  rises  from  9%  (at  2  m  and  20  m)  to 
25%  at  200  m  if  we  assume  no  deposition  between  2  m  and  20  m.  Integrated 
momentum  rises  from  1.3  x  10^  dyne-sec  to  1.1  x  10^  dyne-sec.  While  we 
should  not  put  too  much  stock  in  the  exact  numbers  from  this  model,  it 
seems  clear  that  distant  fireball  coupling  may  not  be  negligible  and  should 
be  studied  with  some  care.  In  section  V  below  we  compare  fireball  coupling 
to  the  effects  of  overpressure  of  the  heated  air,  an  effect  which  may 
dominate  all  the  effects  of  energy  radiatively  coupled  into  the  ground 
outside  of  the  directly  illuminated  region. 

4. 4  Competition  Between  Hydrodynamics  and  Radiation  Diffusion 

In  assuming  constant  p  in  the  equation  for  opacity  we  may  have 
underestimated  the  radiation  coupling,  especially  at  large  radii,  so  we 
must  now  return  to  the  blowoff  question. 

The  rarefaction  wave  of  blowoff  propagates  into  the  ground  with 
the  (hot)  sound  speed  vs  ■  (CyTg)^^  ,  the  time  tjj  at  which  it  catches 
up  with  the  diffusion  wave  is  roughly 

Vh  '  (tHk)1/2  <4-n) 

which  gives 
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We  must  also  compute  t2  ,  defined  as  the  time  that  it  takes  to  heat  a 
single  mean  free  path  X  to  temperature  Tg  »  i«e.,  the  time  at  which  the 
radiation  diffusion  approximation  first  starts  to  be  valid: 


(t2k) 


-  X 


(4.13) 


which  gives 


KT3ac 


(4.14) 


Additionally,  we  need  the  condition  for  a  rarefaction  wave  at  some 
temperature  Tj  <  Ts  to  propagate  through  one  mean  free  path  in  time  tj 
before  photon  heating  exceeds  the  temperature  ,  i.e., 


ApCvTi 


acT4t . 
s  1 


(4.13a) 


and 


X 


<CvV 


(4. 15b) 


Eliminating  Tj  between  these  equations,  we  get 


t 


1 


(4.16) 


Now  eliminating  tj  between  Eqs.  (4.15a)  and  (4.15b)  and  requiring 
Ti  <  Tg  ,  we  obtain 


1 


I 


t5/2  <  p_  c3/2  „  (15(J  eV)5/2  . 


(4. 17) 


What  is  the  meaning  of  all  these  conditions?  In  the  T  plane 
they  map  out  the  different  regimes  of  competition  between  photon  streaming, 
diffusion,  and  blowoff.  Figure  4.4  illustrates  this.  The  various  boundary 
lines  all  meet  (approximately)  at  a  critical  point 
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(4.18) 


In  region  (B)  of  the  figure  diffusion  goes  through  many  mean  free 
paths,  then  the  source  turns  off,  and  only  subsequently  does  the 
rarefaction  wave  catch  up  to  the  diffusion  front. 


In  region  (C)  diffusion  goes  through  many  mean  free  paths,  but  the 
rarefaction  catches  up  while  the  source  is  still  on.  One  then  goes  over  to 
a  regime  in  which  a  coupled  rarefaction-diffusion  wave  penetrates  into  the 


itter  faster  than  t*/2  >  but  slower  than  t  . 


In  region  (D),  matter  starts  to  blow  off  even  before  it  reaches 
the  temperature  Tg  ,  but  it  becomes  heated  to  Tg  during  the  subsequent 
blowoff.  In  this  region  a  coupled  rarefaction-diffusion  similarity 
solution  is  probably  a  reasonable  approximation  during  the  whole  time  the 
source  is  on. 
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Figure  4.4  THE  "TRAJECTORY"  OF  THE  RADIATION  -  GROUND  COUPLING  PROBLEM 
IN  TEMPERATURE-TIME  SPACE.  IN  REGION  (B)  RADIATION  DIFFUSES  MANY  MEAN  FREE 
PATHS  BEFORE  SOURCE  TURNOFF,  HYDRODYNAMICS  IS  NEGLIGIBLE.  IN  (C)  RADIA¬ 
TION  DIFFUSES  MANY  MEAN  FREE  PATHS  BUT  HYDRODYNAMIC  MOTION  CATCHES  UP 
BFFORE  SOURCE  TURNOFF.  IN  ID)  HYDRODYNAMIC  MOTION  OCCURS  BEFORE  GROUND 
MATERIAL  HAS  BEEN  HEATED  TO  THE  SOURCE  TEMPERATURE. 


I 
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Regions  (A)  and  (E)  are  never  of  interest  in  our  problem,  since 
their  times  are  too  small.  In  these  regions,  the  source  time  is  so  short 
that  it  never  heats  even  a  single  mean  free  path  to  its  temperature.  In 
(E),  the  top  mean  free  path  blows  off  while  the  source  is  still  on;  in  (A) 
the  source  turns  off  first. 


Also  in  Figure  4.4,  we  have  indicated  the  locus  of  the  curve 
T(r),t(r)  ,  parametrized  along  its  length  by  r  in  meters.  One  sees  that 
the  curve  enters  region  (D)  at  r  =  35  m  .  Thus  the  distant-fireball 
coupling,  which  we  have  already  noted  as  possibly  important,  may  also  be 
enhanced  by  the  decrease  of  optical  depth  with  rarefaction. 


We  can  estimate  very  roughly  the  possible  enhancement  as 
follows:  if  the  penetration  was  entirely  diffusive,  it  would  go  to  a  depth 


£diff  ~  (tk) 


1/2 


(4.19) 


If  it  were  entirely  hydrodynamical,  it  would  go  to  a  depth 


‘hydro  ~  <W‘/2t 


(4.20) 


The  ratio  ^  £diff  can  seen  ^roni  Eq*  4.11)  to  be  about 

(t/tH)l/2  .  Since  the  true  case  should  be  a  similarity  solution  in  between 
these  two  extremes,  the  enhancement  over  the  predictions  in  Fig.  4.3  should 
be  (t/tH)a^2  where  0  <  a  <  1  .  A  plausible  guess  is  a  «  1/3  to  1/2  , 
which  (from  Figure  4.4)  gives  enhancement  factors  of  10  to  100  for  energy 
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and  momentum  deposition  at  radii  of  50  m.  Clearly  this  is  too  large,  since 
it  would  push  radiation  coupling  efficiencies  to  order  unity.  We  must 
conclude,  however,  that  a  detailed  calculation,  including  details  of  the 
blowoff,  the  change  in  opacity  with  blowoff,  and  radiation  transport  both 
in  the  blown-off  material  and  in  the  structure  of  the  fireball,  is 
necessary. 
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V  HYDRODYNAMICS 

The  purpose  of  this  section  is  to  discuss  the  hydrodynamic 
consequences  of  fireball-coupled  energy  and  momentum,  and  to  describe  a 
simple  model  for  the  eventual  formation  of  a  crater  from  the  energy  and 
momentum  coupled  into  the  ground.  These  estimates  cannot  supplant  the 
results  of  code  calculations;  rather,  their  purpose  is  to  identify  and 
explain  the  dominant  physical  processes  of  these  code  calculations. 

5.1  Blowoff 


X 

3 


% 
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5.1.1  Momentum  Amplification 

A  surface  nuclear  burst  directly  heats  a  very  thin  layer  on  the 
surface  of  the  ground.  Close  in  to  the  burst  the  direct  radiation  heats  to 
a  depth  of  tens  of  centimeters;  in  the  fireball  regime  we  see  from  Fig.  4.2 
of  the  preceding  section  that  the  depth  heated  is  much  less  (an  accurate 
number  requires  a  detailed  calculation).  These  surface  layers  are  raised 
to  very  high  pressure  when  they  are  heated,  and  the  recoil  of  their  blowoff 
delivers  a  sharp,  strong  impulse  to  the  ground  which  sends  a  shock  wave 
propagating  downward. 

The  one  dimensional  problem  of  the  propagation  of  a  strong  shock 
from  a  free  surface  into  a  half  space  is  dealt  with  by  Zeldovich  and 
Raizer^S.  Similarity  solutions  for  these  shocks  in  polytropic  gases  exist, 
and  take  a  simple  form  for  y  ■  7/5  .  This  value  of  y  is  a  fair 
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approximation  to  the  properties  of  ionized  gases  over  a  wide  range  of 
temperature  and  density,  down  to  shock  pressures  of  about  1  Mbar.  For 
y  =  7/5  these  authors  find  the  shock  speed  u  «  t”^3  ,  its  distance 


travelled 

3/5 

x  «  t  ,  its  pressure 

p  .  t"4/5  a  x*4/3  ,  and 

its  impulse  per 

unit  area 

I  *  t1/5  «  X1/3  a  p-1/4 

;  other  relations  may  be 

trivially 

derived  (only  one  relation  is  independent).  What  is  striking  about  these 
results  is  the  increase  in  delivered  impulse  with  the  depth  of  penetration 
of  the  shock.  This  phenomenon  has  been  colloquially  called  "momentum 
amplification."  As  the  shock  penetrates  into  the  ground  more  and  more 
matter  blows  off,  carrying  more  and  more  momentum  upward,  and  producing  a 
correspondingly  increasing  downward  recoil.  Equivalently,  the  energy  in 
the  ground  is  shared  by  an  increasing  amount  of  matter  necessarily  at  a 
decreasing  velocity.  The  increasing  amount  of  matter  involved  means  that  a 
given  amount  of  energy  corresponds  to  an  increasing  impulse  in  the  ground 
(balanced,  of  course,  by  the  momentum  carried  by  the  blown  off  material). 
This  phenomenon  is  complicated  by  the  continuing  loss  of  energy  from  the 
ground  to  the  blowoff,  a  loss  whose  importance  depends  on  the  equation  of 
state.  For  a  soft  (isothermal)  equation  of  state  y  *  1  and 
I  «  t°  ;  there  is  no  momentum  amplification  because  the  energy  in  the 
downward  moving  ground  and  its  velocity  decrease  at  the  same  rate.  For  the 
hardest  polytropic  equation  of  state  y  ♦  “>  and  1  «  j.0.284  ^  slightly 
more  momentum  amplification  than  for  y  ■  7/5  because  now  the  energy  in 

downward  moving  ground  decreases  very  slowly  («  t”’^^)  while  the 

—  358 

velocity  decrease  u«  t  *  is  still  rapid. 


Close  in  to  the  burst  the  deposition  of  radiative  energy  is 
rapidly  varying  horizontally,  so  the  ground  motion  is  a  two  dimensional 
problem.  Fortunately,  this  problem  is  well  calculated  by  existing  codes, 
and  is  not  controversial.  There  is  more  argument  about  the  fate  of  the 
coupled  energy  and  momentim  in  the  fireball  region,  where  the  motion  is 
very  close  to  slab-symmetric,  but  where  the  initial  deposition  is  in  such  a 
thin  layer  that  numerical  calculation  may  be  tricky.  Although  millimeters 
of  solid  ground  are  eventually  heated,  the  first  mean  free  path,  which  is 
the  thickness  of  the  initial  blowoff  in  most  of  the  fireball,  is  only  about 
lOvm  thick  (see  Fig.  4.2).  It  is  here  that  we  wish  to  discuss  the 
momentum  amplification  problem  using  the  results  of  Zeldovich  and  Raizer. 

The  inner  parts  of  the  fireball  (r  <  35  m)  are  in  region  C  of 
Fig.  4.4  of  the  preceding  section.  Initially  the  fireball  heats  at  least 
one  optical  depth  to  its  effective  temperature  T  >  150  eV  .  At  solid 
density  this  hot  ground  has  a  pressure  ~  100  Mbar  ,  and  the  equation  of 
state  is  that  of  a  gas.  Under  adiabatic  conditions  it  may  be  reasonably 
described  with  a  polytropic  equation  of  state  P  «  ,  where  y  lies 

between  1  and  5/3,  depending  on  how  rapidly  the  degree  of  ionization  is 
varying  with  temperature,  and  hence  on  entropy.  We  shall  use  y  =  7/5 
throughout. 

Even  for  the  most  favorable  equations  of  state  momentum 
amplification  is  a  slow  process.  It  is  limited  by  the  range  in  pressure 
over  which  the  matter  retains  its  polytropic  equation  of  state.  For  shock 
pressures  below  1  Mbar  the  matter  is  not  in  a  vapor  or  plasma  phase  upon 


unloading,  but  rather  is  a  fractured  solid,  the  hydrodynamics  of  which  are 
much  closer  to  that  of  an  incompressible  fluid  than  that  of  a  polytropic 
gas.  A  pressure  of  1  Mbar  is  the  lower  bound  on  the  region  of  validity  of 
the  momentum-amplifying  blowoff,  so  the  total  amplification  for  y  =  7/5 
is  (100)^^  •  3  .  In  practice  momentum  amplification  is  even  smaller  than 
this  for  two  reasons:  (1)  porosity  and  microstructure  (sand,  snow, 
vegetation,  etc.)  in  the  uppermost  layers  of  ground  reduce  the  density,  and 
hence  the  pressure,  at  which  the  amplifying  process  begins  and  (2)  no 
momentum  amplification  applies  to  energy  delivered  after  the  initial 
rarefaction  catches  up  to  the  radiation  diffusion  front,  because  such 
energy  is  not  deposited  in  solid-density,  high  pressure  (100  Mbar)  matter. 

Beyond  35  meters  (standard  fireball),  the  fireball  radiation  is 
softer  (T  <  150  eV),  and,  as  discussed  in  the  preceding  section, 
penetrates  less  efficiently  into  the  ground.  Most  of  the  fireball  coupling 
occurs  in  the  range  between  50  and  100  meters,  where  T  <  100  eV  .  No 
portion  of  the  ground  is  ever  heated  to  the  fireball  temperature  before  it 
is  cooled  by  rarefaction.  Peak  pressures  are  then  much  lower  than  regions 
for  which  T  >  150  eV  ,  partly  because  the  fireball  temperature  is  lower 
and  partly  because  the  ground  does  not  reach  even  this  reduced  temperature. 
In  such  regions  there  is  reduced  opportunity  for  momentum  amplification. 

The  fundamental  reason  why  momentum  amplification  over  most  of  the 
fireball  region  is  insignificant  may  be  seen  from  the  time  scales  in  Fig. 
4.4.  For  T  <  150  eV  the  fireball  duration  exceeds  tj  ,  the  time  for  the 
rarefaction  to  catch  up  with  the  heating  front,  by  several  orders  of 
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magnitude.  Only  a  lOpm  thick  skin  is  heated  during  the  first  interval 
.  For  the  case  T  >  150  eV  which  applies  only  over  a  small  region  the 
fireball  duration  should  be  compared  to  tjj  ;  the  inequality  is  still 
large.  Energy  deposited  after  the  brief  initial  time  tj  f  (or  tH),  leads 
to  little  or  no  momentum  amplification,  because  it  is  deposited  in  low 
density  blowoff  at  pressures  far  below  the  peak  of  100  Mbar.  This  later 
time  includes  more  than  99.9%  of  the  fireball  deposition  outside  30  m,  the 
beginning  of  the  fireball  regime.  Even  the  existence  of  the  initial 
transient  is  dubious,  because  the  front  of  a  real  fireball  does  not  produce 
an  instantaneous  jump  in  its  radiation  field,  but  has  some  finite  width. 

5.1.2  Air  Tamping 

After  the  rarefaction  front  catches  up  to  the  thin,  initially 
fireball-heated,  surface  layer,  a  coupled  rarefaction-thermal  front 
penetrates  into  and  ablates  the  ground.  As  discussed  in  the  preceding 
section,  it  is  hard  to  calculate  simply  with  pencil  and  paper  the  progress 
of  the  front  and  its  depth  of  penetration,  but  it  is  tractable  on  a 
computer.  Some  general  considerations  are  easy  to  establish,  however.  The 
ablated  matter  is  heated  to  the  fireball  temperature,  but  only  after  it 
rarefies.  The  maximum  pressure  is  far  below  the  100  Mbar  value  which  solid 
density  matter  has  at  fireball  temperature,  so  that  the  maximum  possible 
momentum  amplification  is  reduced  below  the  bound  of  3;  it  cannot  be  much 
more  than  unity.  If  the  blowoff  were  into  vacuum  the  maximum  density  (and 
hence  maximum  pressure)  in  the  flow  would  be  a  steadily  decreasing  function 
of  time,  as  a  thickening  layer  of  ablated  material  built  up  between  the 
radiation  source  and  the  solid  ground.  Reality  is  different  (until  we 


fight  wars  on  the  Moon).  The  ground  is  covered  with  air  at  fireball 

temperature.  At  a  density  of  10-3  g/cm3  and  at  fireball  temperatures,  air 

is  multiply  ionized  (generally  five  times  for  nitrogen  and  six  times  for 

oxygen)  and  has  a  pressure  of  ~  10-30  Kbar.  Once  the  ground  blowoff  has 

rarefied  to  the  same  density  as  air  the  pressures  of  the  two  materials  are 

nearly  equal  (both  are  multiply  ionized,  and  have  similar  mean  Z  ).  Ibis 

happens  rather  quickly  in  a  relaxation  time  tr  of  about  p  /p  .  times 

t  g  air 

the  initial  rarefaction  time  tj  or  tH  •  At  the  100  eV  point  in  our 
standard  fireball,  t  ~  10~^  which  corresponds  to  about  10~3  Qf  the 
fireball  duration.  Then  the  air  decelerates  the  expanding  front  of  the 
blowoff.  After  a  time  equal  to  a  few  times  tr  ,  pressure  pulses  are 
exchanged  through  the  blowoff  from  air  to  ground  and  back  again,  and  as 
transient  effects  die  away  the  blowoff  air  and  upper  portion  of  the  dense 
ground  settle  down  to  a  nearly  isobaric  condition.  The  reason  for  this  is 
the  short  relaxation  time  for  pressure  through  the  thin  layers  between  the 
semi-infinite  atmosphere  (already  isobaric)  and  the  semi-infinite  dense 
ground.  After  these  initial  transients,  all  other  properties  of  the 
ablating  flow  change  slowly.  Eventually  (in  lCT^  -  10-3  s)  this  steady 
state  is  brought  to  an  end  by  the  turning  off  of  the  fireball,  or  the 
progress  inward  of  a  rarefaction  in  the  air  from  its  upper  surface.  Until 
then,  the  fireball  and  ablating  surface  act  as  a  simple  pressure  boundary 
condition  on  the  ground,  and  a  planar  shock  propagates  downward.  Because 
this  shock  is  weak  (the  overpressure  is  tens  of  kilobars)  it  produces 
negligible  compression,  no  blowoff,  and  no  vaporization. 


It  may  be  seen  from  Fig.  5.1  that  the  "airslap"  contribution  to 
ground  impulse  we  have  just  described  appears  to  be  dominant  (remember, 
though,  that  for  the  close  in  coupling  momentum  amplification  may  be  a 
factor  of  ~  10).  Particularly  striking  is  the  fact  that  the  delivered 
impulse  is  still  rapidly  increasing  at  200  meters,  the  limit  of  available 
fireball  parameters.  What  limits  this  impulse? 

Airslap  is  large  because  area  and  dwell  time  both  increase  rapidly 
outward  in  the  fireball,  while  temperature  decreases  only  moderately. 
Eventually  though,  and  not  too  far  beyond  200  m,  the  fireball  reaches  its 
maximum  extent;  beyond  this  point  its  dwell  time  drops  abruptly  to  zero, 
and  there  is  no  fireball.  There  is  still  some  airslap  from  shocked  air, 
but  its  overpressure  is  negligible.  A  second  limiting  factor  enters  at  the 
outer  edge  of  the  fireball.  The  air  overpressure,  which  is  about  10  Kbar 
at  200  meters,  is  dropping.  Once  it  drops  below  the  crushing  strength  of 
rock  (one  or  a  few  Kbar)  the  shocked  rock  can  unload  elastically,  and  the 
delivered  impulse  does  not  contribute  to  cratering.  Because  of  the  large 
contribution  to  impulse  coupling  of  the  outer  edge  of  the  fireball,  it  is 

important  that  computer  codes  treat  the  late  time,  low  pressure  regime  well 

if  accurate  results  are  to  be  obtained. 

5.2  Incompressible  Hydrodynamics 

Under  the  fireball  is  a  weak,  downward  propagating  shock.  Its 

peak  pressure  is  much  less  than  1  Mbar,  so  it  produces  very  little 

compression.  It  travels  at  a  velocity  very  close  to  the  compressional 

sound  speed  v  (v  =  few  km/s)  ;  the  shocked  material  has  a  subsonic 
s  s 
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Figure  5.1  CONTRIBUTIONS  TO  THE  MOMENTUM  IMPULSE  GIVEN  TO  THE  GROUND 
FROM  THE  RADIATION  INDUCED  VACUUM  BLOWOFF  DESCRIBED  IN  THE 
PRECEDING  SECTION  AND  THE  AIRSLAP  DISCUSSED  IN  THIS  SECTION 
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downward  velocity  of  v  *>  — —  which  is  ~  100-300  m/s  .  These  conditions 

pVs 

continue,  the  shock  propagating  downward,  until  the  fireball  ceases  and  the 
surface  pressure  relaxes  to  zero.  Then  a  rarefaction  follows  the  shock 
leaving  in  its  wake  matter  at  rest,  but  displaced  downward.  In  slab 
symmetry  the  rarefaction  would  chase  the  shock  deep  into  the  earth,  very 
slowly  attenuating  it.  In  the  real  world  the  fireball  is  of  finite 
hori.  .ital  extent.  The  shock  propagates  sideways  in  addition  to  downward, 
and  as  it  propagates  horizontally  the  free  surface  is  a  continuous  source 
of  rarefactions.  When  rarefactions  interact  they  produce  a  state  of  net 
tension  in  which  the  rock  fragments.  The  final  displacement  of  this  loose 
material  creates  the  crater. 

If  the  eventual  crater  is  larger  than  the  fireball,  as  we  expect, 
its  size  may  be  roughly  estimated.  All  of  the  momentum  delivered  to  the 
ground,  whether  close  in  to  the  burst  or  out  at  the  limits  of  the  fireball, 
may  be  combined  into  a  single  downward  impulse  concentrated  at  one  point. 
The  details  of  compressible  hydrodynamics  near  the  burst  and  at  the  surface 
under  the  fireball  (pressures  over  1  Mbar,  or  temperatures  over  1  eV)  are 
ignored — all  that  matters  is  the  impulse  they  deliver  into  a  much  larger 
volume  which  undergoes  essentially  incompressible  flow.  Because  the 
stresses  far  exceed  the  yield  strength  of  rocks  we  may  assume  that  the  flow 
is  hydrodynamic,  until  the  outer  limits  of  the  cratered  volume  are  reached. 
As  sketched  in  Fig.  5.2  the  downward  impulse  I  is  spread  over  a  steadily 
widening  area.  Approximately  we  can  write 

2 

I  ■  PR  At  (5.1) 
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where  P  is  a  typical  pressure  in  the  flow  and  At  its  duration.  Then, 
similarly  approximately,  At  =  R/v  ,  where  v  is  a  typical  velocity, 
and  M  is  the  mass  swept  up  in  the  flow 


v  =  I/M  =  I/pR3  .  (5.2) 

The  size  of  the  cratering  flow  is  thus  estimated  by 

R  =  I1/3(pP)"1/6  .  (5.3) 

As  the  flow  spreads  out,  and  the  momentum  is  spread  over  an 
increasing  mass  of  crushed  and  fragmented  rock,  kinetic  energy  is  lost  to 
internal  and  external  friction  and  to  fracture  energy.  This  process 
terminates  when  P  reaches  the  strength  of  the  rock;  lower  stresses  are 
carried  away  by  elastic  waves.  Then  the  final  crater  size  is  approximately 
given  by  R  «  I^^pP)-*^  .  beyond  this  radius  the  rock  is  unfractured. 
This  estimate  is  crude,  and  is  particularly  bad  for  unconsolidated 
materials  like  sand  whose  strength  is  strongly  dependent  on  loading.  As 
before,  our  estimates  are  not  meant  as  a  substitute  for  computer 
calculations,  but  rather  as  a  simple  model  to  assist  understanding. 
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VI  CONCLUSIONS  AND  RECOMMENDATIONS 

(1)  We  have  found  the  cratering  problem  to  be  a  complex  one 
indeed.  Its  intrinsic  complexity  is  exacerbated  by  the  poorly  documented 
experimental  data  for  the  few  existing  surface  bursts  of  modern,  radiating 
nuclear  weapons  and  by  the  small  overlap  in  source  configurations,  material 
descriptions,  etc.,  which  have  been  employed  in  computer  studies.  None  of 
this  is  new  to  the  workers  in  the  field  who  indeed  seem  to  share  the  same 
frustrations.  Further  efforts  to  compare,  contrast  and  refine  theoretical 
calculations  are  strongly  encouraged. 

(2)  The  apparent  success  of  computer  predictions  for  crater 
formation  by  high  explosives  (HE)  in  simple  geologies  is  highly  relevant 
since  one  must  work  through  similar  hydrodynamic,  elasto-plastic  and  rock 
fracture  processes  in  the  nuclear  case  as  well.  If  the  HE  success  is  real 
then  it  follows  that  if  one  can  calculate  accurately  the  energy  and 
momentum  configuration  in  the  ground  at  the  end  of  the  radiation  era,  one 
ought  to  be  able  to  carry  the  remainder  of  the  calculation  through  to  a 
reliable  result  for  a  nuclear  crater.  Comparisons  of  code  and  centrifuge 
HE  cratering  experiments  have  not  yet  been  done  but  would  provide  valuable 
checks  on  numerical  calculations  in  well  defined  situations. 

(3)  The  crude  model  of  radiation  coupling  we  develop  in  section 
IV  gives  order  of  magnitude  confirmation  of  the  predictions^?*^  that  only 
a  small  fraction  (few  percent)  of  the  direct  radiation  is  coupled  into  the 
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ground.  This  directly  coupled  energy  is  concentrated  in  a  small  region 
(r  <  5  m)  and  occurs  during  t  <  100  ns  for  a  1  MT  burst.  The  predicted 
depth  to  which  this  radiation  penetrates  is  also  qualitatively  reproduced 
by  our  model  (d  =  60  cm)  . 

(A)  The  radiation  model  of  section  IV  also  indicates  that  the 
total  (integrated  over  area)  energy  which  is  coupled  from  the  fireball  into 
the  ground  at  large  distances  from  the  source  (r  «  40-200  m),  is  comparable 
to  the  directly  coupled  energy.  This  comes  about  in  spite  of  the  lower 
fireball  temperature  because  of  the  long  duration  of  the  fireball  at  these 
radii  (t  »  10  -  10  s)  and  the  quadratic  increase  of  area  with 

radius.  The  initial  impulse  generated  by  the  fireball  coupling  is  also 
predicted  to  be  comparable  to  the  initial  impulse  generated  in  the  close-in 
region  by  the  directly  coupled  radiation.  The  depth  of  the  fireball 
heating  is,  however,  very  small.  Precise  values  depend  on  the  solution  of 
a  combined  hydrodynamic  blowof f-radiation  diffusion  problem  we  have  not 
undertaken.  At  r  ■  50  m  the  heating  depth  is  expected  to  lie  somewhere 
between  0.4-4  cm  and  to  decrease  to  0.01-0.1  cm  at  r  =  200  m. 

(5)  We  are  not  sure  if  the  net  effect  of  fireball  coupling  will 
be  sizable  enough  to  have  a  significant  influence  on  final  crater  size  or 
configuration.  We  believe,  however,  just  as  for  the  problem  of  direct 
radiation  coupling,  that  the  correct  answer  is  fully  within  the  grasp  of 
present  theoretical  knowledge  and  computer  capability. 

(6)  In  an  effort  to  effect  a  convergence  among  the  various 
numerical  code  results  for  the  radiation  coupling,  we  recommmend  a  series 
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of  bench  mark  calculations.  These  might  be  designed  as  follows:  In  order 
to  separate  the  complications  of  geometry  from  the  rad iation/hydrodynmamic 
issues  of  primary  interest,  a  series  of  one  dimensional  calculations  could 
be  undertaken  which  correspond  to  a  wide  variety  of  radiation  temperatures 
(or  spectral  conditions  as  appropriate)  and  time  durations.  The  range  of 
these  parameters  should  cover  conditions  from  those  appropriate  to  the 
direct  radiation  all  the  way  to  those  appropriate  to  typical  fireball 
conditions  at  late  times.  In  this  way,  one  could  develop  a  clear  sense  of 
the  various  depth  and  time  scales  involved  in  all  parts  of  the  cratering 
problem  at  all  radii.  This  knowledge  could  then  be  used  to  assure  that 
nothing  important  is  being  lost  in  the  full  geometry,  realistic  crater 
calculations. 

(7)  As  part  of  the  radiation  studies  we  recommend  that 
considerable  attention  be  given  to  the  sensitivity  to  height  of  burst. 
Since  it  may  very  well  be  that  modern  nuclear  weapons  do  indeed  produce 

f  small  craters,  the  dependence  on  height  of  burst  may  be  complex. 

Calculations  over  a  wide  range  of  heights  are  needed.  In  doing  this  work, 

/ 

V  the  source  should  be  characterized  in  simple  terms  rather  than  including 

I  complex  details  of  design  and  configuration. 

(8)  When  satisfactory  convergence  has  been  obtained  on  the 

^  radiation  coupling  calculations,  and  especially  if  the  results  settle  on 

small  efficiencies  of  cratering,  one  ought  then  to  go  back  and  trace 
through  the  fate  of  the  other  (small)  weapon  energy  releases — case  debris, 
I  neutron  emission,  etc. 
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(9)  It  may  be  particularly  important  to  accurately  calculate  the 
effects  of  airslap  in  the  outer  parts  of  the  fireball,  and  to  understand 
the  yield  properties  of  the  undisturbed  ground  at  stresses  not  far  above 
its  elastic  limit. 

(10)  Since  it  seems  to  be  clear  that  the  mechanisms  of  coupling 
and  the  resulting  spatial  distributions  for  modern  nuclear  weapons  are 
quite  different  from  those  of  high  explosive  bursts,  the  same  scaling  laws 
need  not  apply  to  both.  It  could  therefore  be  very  misleading  to  apply  the 
empirical  scaling  laws  which  have  been  extracted  from  the  extensive  HE  data 
base  to  nuclear  weapon  bursts. 

(11)  There  is  evidently  a  considerable  disagreement  over  the 
appropriate  description  of  material  properties  and  yield  criteria  which  go 
into  cratering  calculations  and  the  level  of  detail  which  is  necessary.  It 
would  be  valuable  to  have  some  simplified  reference  calculations  which 
would  indicate  the  sensitivity  to  various  descriptions. 

(12)  We  have  little  to  say  about  the  craters  which  were  produced 
during  testing  in  the  Pacific  Proving  Ground  (PPG).  The  actual  device 
environments  are  so  unknown  and  the  crater  modifications  (wave  washing, 
slumping,  etc.)  which  may  have  occurred  after  the  burst  are  so  uncertain, 
that  the  PPG  nuclear  data  may  have  been  rendered  essentially  useless. 
Applications  to  the  PPG  craters  of  scaling  laws,  which  themselves  have  very 
limited  verification,  should  be  treated  with  great  caution. 
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